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Abstract. Power suppression of the cosmic microwave background on the largest observable 
scales could provide valuable clues about the particle physics underlying inflation. Here we 
consider the prospect of power suppression in the context of the multifield landscape. Based 
on the assumption that our observable universe emerges from a tunnelling event and that 
the relevant features originate purely from inflationary dynamics, we find that the power 
spectrum not only contains information on single-field dynamics, but also places strong con¬ 
straints on all scalar fields present in the theory. We find that the simplest single-field models 
giving rise to power suppression do not generalise to multifield models in a straightforward 
way, as the resulting superhorizon evolution of the curvature perturbation tends to erase 
any power suppression present at horizon crossing. On the other hand, multifield effects do 
present a means of generating power suppression which to our knowledge has so far not been 
considered. We propose a mechanism to illustrate this, which we dub flume inflation. 
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1 Introduction 

There is substantial evidence both on theoretical and observational fronts indicating that 
our observable universe underwent a period of inflation [1-4]. If inflation really occurred, it 
provides an extraordinary opportunity to use observations of the Cosmic Microwave Back¬ 
ground (CMB) and large scale structure as a means of studying particle physics at energy 
scales well above those likely to be achieved with terrestrial experiments. In particular, in¬ 
flation is sensitive to Planck suppressed operators, allowing us to test ultraviolet complete 
theories in a unique way. It remains a central challenge for inflation to understand how it 
is to be embedded in such a theory but there has been significant progress in the context of 
string theory (for a recent review see [5]). 

Developments in the understanding of compactification over the last decade have given 
rise to a striking picture that should have profound implications for the study of inflation — 
the prevalence of large Hodge numbers [6-9] indicates the existence of many, often hundreds of 
scalar fields, interacting via a complicated potential containing a large number of metastable 
vacua. This picture is sometimes referred to as the string landscape [10, 11]. Eor inflationary 
phenomenology, arguably the most obvious consequences of this scenario are the possibility 
of multifield dynamics and the idea that our observable universe might have originated from a 
tunnelling event. Both of these behaviours present the alluring prospect of observable effects 
in the CMB. 
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The tunnelling process gives rise to an open universe inside the bubble of the new 
vacuum, whose initial evolution is dominated by spatial curvature [12], This leaves an imprint 
on the power spectrum at scales crossing the horizon during this period [13-17] If the 
inflationary epoch inside the bubble is short enough, one could hope to see this effect in the 
low-l' (largest observables scales) regime of the CMB spectrum. This possibility is heavily 
constrained by the latest Planck data [23], which sets a bound on the scale of curvature as 
being at least an order of magnitude bigger than the scale of today’s horizon, implying that 
the scales we observe left the horizon after curvature domination. 

Multiheld effects are potentially very important for inflation embedded in fundamental 
physics. Obtaining an extended period of inflation is notoriously difficult in string theory 
but it is often the case that whatever mechanism enables one field to be sufficiently light, will 
tend to make other fields sufficiently light to be cosmologically relevant as well. This point 
has been emphasised in a large body of work but the message in Ref. [24] is particularly 
succinct; multifield dynamics should be considered the norm in string theory. 

Inflation with more than one light field can in principle give rise to a wide range of 
signatures, generally (but not exclusively) as a consequence of the superhorizon evolution 
of the primordial curvature perturbation, which in turn is a consequence of the presence 
of isocurvature degrees of freedom. Currently there is no observational evidence for any of 
these signatures, implying no evidence of multifield inflation. However it should be noted 
that this does not provide a confirmation against these scenarios either, since many models 
that exhibit multifield dynamics do not break the single-field predictions enough to be under 
threat from current constraints (see for instance [25-33]). 

In this work we will focus on the particularly interesting observation of a possible 
suppression of the power spectrum on large scales, as suggested by the WMAP [34] and 
Planck [35] data. There has been considerable interest in studying mechanisms for this lack 
of power [36-52], in particular as a consequence of the inflationary era, but its possible mi¬ 
crophysical origin remains undisclosed. It was pointed out in Ref. [53] that tunnelling from 
a metastable vacuum through a potential barrier might naturally generate a period of steep 
inflaton potential before the slow-roll plateau. Should we be fortunate enough that the total 
period of inflation is sufficiently small and, at the same time, long enough to avoid the con¬ 
straints of curvature, we could hope to see the onset of this epoch imprinted in the CMB as 
some form of signature on large scales. Of key relevance to this paper, Bousso, Harlow and 
Senatore [40, 45] emphasised that in the case of a theory consisting of a single scalar field, 
power suppression is a likely consequence of this tunnelling in the landscape. 

The possibility that power suppression may be a consequence of the string landscape is 
extremely exciting, as it would represent a unique door to the physics of the early universe. 
However there is a considerable way to go before such claims can be made. Building on the 
work of Ref. [40, 45], an obvious next question to ask in furthering this pursuit is does power 
suppression also occur if more than one field is involved in this event? 

The tunnelling process in a multifield model can be quite complicated and one might 
expect the initial condition after the transition to be quite far from the inflationary plateau 
region that leads to most of the observable scales in the CMB, the ^ > 50 range One 

^Another possible observational consequence of this scenario is the detection of collisions between our 
bubble and other topological defects [18] or bubbles [19-22] created in the parent vacuum. This is an interesting 
possibility that is currently being actively investigated. 

^See [54] for a proposal that links the position of the fields after the tunnelling event to the inflationary 
region in a multifield landscape. 


- 2 - 



can expect scenarios where many moduli fields are affected by the tunnelling transition, 
potentially giving rise to a rich mass spectrum of the scalar fields involved in inflation. As 
an example of hierarchies in the mass spectrum, one can conceive that tunnelling would 
affect different moduli sectors in distinct ways. This could result, for example, in a situation 
where at first inflation is controlled by the evolution of complex structure moduli and then, 
as these fields settle into their minima, becomes dominated by a different moduli sector, 
like a Kahler moduli field. Many current models of compactification naturally lead to this 
kind of hierarchies in the mass spectrum and separation in scalar sectors, indicating that 
multifield dynamics are relevant after tunnelling events (see, for example Ref. [6]). In this 
paper we explore how the single-field results of power suppression on large scales generalise 
in the presence of many scalar helds with such rich mass hierarchies. 

We start our discussion by reviewing a simple method to compute the power spectrum 
and the spectral index in a multifield inflationary context. We then review the mechanism of 
power suppression due to a steepening of a single-field potential as described in Ref. [40, 45], 
and look at the consequences of including more than one light field in this scenario. Our 
main result is that achieving power suppression is considerably more delicate in a multifield 
model than in the single-field case. In fact, a steep potential and a fast evolution during the 
first few e-folds of inflation do not necessarily imply suppression of the power spectrum on 
large scales. We identify strong constraints imposed on initial conditions as well as inflaton 
potential in the region resulting in suppression. These constraints apply not only to the form 
of the potential along the inflationary trajectory but to all other directions in field space. 
This is a rare opportunity where 2-point statistics can be used to learn about multifield 
effects. 

If all relevant features originate purely from an inflationary era, suppression on large 
scales would then constrain multifield dynamics in a region where one also requires non-trivial 
evolution of the slow-roll parameter e. We argue that the ability to learn about the mass 
spectrum at such a crucial stage of inflation has significant implications for model building 
in the context of the string landscape. 

To conclude, we present a novel approach to obtain relative power suppression on large 
scales based entirely on superhorizon evolution of the perturbations. Here the idea is not 
to suppress the power at large scales but to enhance it at smaller scales by continuously 
transferring power from isocurvature to adiabatic modes. We present a simple example 
where this can happen which could be relevant for some string theory scenarios. 


2 Multifield infiation — a geometrical picture 

Our goal is to gain intuition about the phenomenology underlying power suppression and so 
to aid us in this task we will compare two methods of computing the perturbations. One 
approach is the transport method — a fully numerical approach which solves for the field- 
field, field-momenta and momenta-momenta correlation functions, which we do not review 
here. We use the publicly available Mathematica code available from transportmethod.com to 
perform our analyses and refer the reader to Ref. [55] for a detailed description of the method. 
Our second approach is a slow-roll superhorizon analysis of the perturbations which we do 
now review. The key benefit of this second method is that it enables one to derive a number 
of analytic or semi-analytic expressions which will be our primary tool in understanding the 
phenomenology underlying our numerical results. 
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Restricting our studies to superhorizon scales and assuming the slow-roll approximations 
radically simplifies the computation of observable quantities. On superhorizon scales the field 
perturbations evolve classically in the sense that decaying solutions to the mode equations 
have died away. The upshot of utilising the slow-roll equations is that the evolution of all 
quantities of interest can be understood purely in terms of gradient flow. Hence a prob¬ 
lem which in principle requires understanding potentially complicated operator equations, 
under these restrictions, reduces to a question of computing relatively simple geometrical 
quantities in field space. This simplification has been used in a very large body of work and 
is intimately related to the celebrated separate universe assumption [56-60]. We refer the 
reader to Ref. [61] and references therein for a more detailed discussion. The notion of using 
geometrical quantities to compute observables was most heavily emphasised in [61-64] and 
it is these works which we summarise here. 

Throughout this paper we work in units of c = li = 1, so that the reduced Planck mass 
reads = SttG = 1, and assume the background space-time metric to be the Friedman- 
Robertson-Walker with “mostly plus” signature (—, -|-, -|-) 

ds^ = —dt^ -|- a{t)‘^dy?. ( 2 . 1 ) 


2.1 Background equations of motion 

We will consider models describing the inflationary dynamics of a set of N real scalar fields, 
i = 1,..., N. Restricting ourselves to theories involving at most two space-time derivatives, 
the action can be written as 

S = - j d^x^g(\G^j dV + , (2.2) 


where Gij = Gij{(j)) is an arbitrary symmetric matrix and V{4>) is the scalar potential. 
The matrix Gij is usually interpreted as a metric on a N— dimensional scalar manifold A4 
parametrized by the fields The equations of motion for the homogeneous background 
fields (/)* = are 


dt 


+ + G^^VjV = 0 , 


(2.3) 


where dots denote derivatives with respect to cosmological time, H = ^ is the Hubble 
parameter, and we have defined the covariant derivative involving the Christoffel symbols 
associated to Gij as DX^ = dX^ + 

In the following we will focus on the evolution of the set of fields that are light enough 
to be part of the slow-roll dynamics. This does not preclude the existence of a heavy sector 
but we will assume that these other fields do not play a significant role for the background 
or the perturbations. However, situations where this assumption fails can be captured using 
the transport method and we do include examples of this in later sections. 

The slow-roll conditions in the multifield case can be defined in terms of the parameters ^ 


e = -H/H^ <C 1; |M| = |Mj| = |V%TogR| <C 1. (2.4) 

Under these conditions, the equations of motion reduce to gradient flow equations 


(/,*' =log U, (2.5) 

where we have used the short hand notation ' = ^ to refer to derivatives with respect to 
the number of e-folds dN = Hdt. 

^We raise and lower indices using the metric Gij and its inverse 
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2.2 Evolution of scalar perturbations 

To study fluctuations about the homogenous background we choose to work in the flat gauge 
such that the independent degrees of freedom are the fluctuations in the fields These 
perturbations transform covariantly under a change of coordinate basis. In other words, 54>^ 
should be interpreted as a tangent vector rather than a coordinate displacement [65]. Our task 
now is to derive an equation of motion for these perturbations. There are at least two ways to 
do this. One method is to perturb the action (2.2) as was done in Refs. [55, 58, 66], however, 
given our simplifying assumptions of only considering superhorizon scales in the slow-roll 
regime, there is a more direct method which is to instead perturb Eq. (2.5). The separate 
universe assumption provides an intuitive picture of why this method works. It states that 
when spatial patches are smoothed on a scale much larger than the horizon scale, the average 
evolution of each patch evolves according to the background equations of motion [56-60]. 
The resulting picture in field space is a bundle of non-interacting trajectories, each evolving 
according to Eq. (2.3) but subject to perturbed initial conditions. When slow-roll holds such 
that Eq. (2.5) applies, this description becomes precisely analogous to geometrical optics [61]. 
Hence, the field fluctuations 6(j)^ may be interpreted as Jacobi fields satisfying [61, 62, 64, 66] 


DScf) 

dN 


-MScj), 


( 2 . 6 ) 


where M can be seen as an effective mass matrix, which encodes the couplings between the 
helds as well as a correction due to the non-trivial geometry of the scalar manifold. Denoting 
^]ki components of the Riemann tensor on A4, the matrix M is dehned as ^ 

M = M - -R with R] = (2.7) 

3 

Ultimately the principal observables of interest are the correlation functions of the primordial 
curvature perturbation We therefore need an expression that relates field perturbations in 
the flat gauge to the primordial curvature perturbation in the constant density gauge ^ . This 
gauge transformation can be computed using cosmological perturbation theory, however once 
again the separate universe assumption enables us to take a shortcut. Lyth and Rodriguez [60] 
showed that the separate universe assumption can be used as a practical means of computing 
Q since on superhorizon scales C = 5N, the variation in the number of e-folds between an 
initial flat slice and a subsequent constant density slice. At lowest order we have 


C(iVf) = 5N 


Hf-. 




N=N{ 


( 2 . 8 ) 


where v = y/Gij((>'^'4>R = and for clarity in discussions to follow, we have explicitly 
labeled the fact that all quantities are to be evaluated at the final time of interest N{. We 

the slow-roll limit we must also require that |f?j | <C 1 
®We choose to work with the curvature perturbation in the constant density gauge, but we could have 
equally worked in the comoving gauge where the curvature perturbation is usually denoted TZ. Both of these 
quantities can be computed using cosmological perturbation theory and are known to be equal at second order 
on superhorizon scales up to 0{k/aH)^ corrections [67, 68]. 
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refer the reader to Ref. [69] for a recent discussion of this topic and detailed derivations of this 
expression as well as higher order expressions using both cosmological perturbation theory 
and the separate universe assumption. 

So far we have implicitly been working on the coordinate basis for our perturbations, 
meaning the basis vectors are aligned with the original fields. One can also look at the vector 
of perturbations Scf) projected onto the so-called kinematic basis [62, 64, 70, 71] defined by a 
set of N orthonormal vectors = {eiue"} satisfying 


eii = 


0' 


~dN 




6a ■ 66 — babj 


(2.9) 


where a = 2,..., N, and we have denoted the dot product between vectors by al • 6 = G^^aibj. 
Note that the quantity Z 21 , which by definition we choose to be non-negative, only vanishes 
when the background follows a geodesic trajectory, i.e. when eyl • Vey = 0. The advantage of 
this basis is that by comparison with Eq. (2.8), one can immediately identify the projection 
of the field perturbations along the inflationary trajectory with the primordial curvature 
perturbation One can then write the following decomposition of the perturbations in the 
kinematic basis corresponding to the time N = Nf 


N 


S(f) = vCeii + y 64>±e' 


a=2 


( 2 . 10 ) 


In this basis, the equations of motion for the superhorizon evolution of the perturbations 
become, ® 

d2) 


= -2Z21^ 


dN 

d6(j)^ 

dN 


( 2 . 11 ) 


where (2) refers to the element a = 
introduced the matrix 


= -[M - Z]„/3 
2, not to be confused with an exponent 2. We have 


■^a6 = 6t 


dN ’ 


( 2 . 12 ) 


that describes how quickly the kinematic basis vectors change along the inflationary trajec¬ 
tory; in other words, it expresses the turn rates of the basis vectors. 

Looking at the equations of motion for the perturbations in this decomposition one 
can easily see two very important points. The first equation tells us that the curvature 

perturbation is not constant, as it can get sourced by isocurvature modes via the first mode 
( 2 ) 

5(jy^ . This will happen whenever there are turns in the trajectory in field space, or in 
other words, deviations from a geodesic motion, as this makes Z 21 non zero. The rest 
of the equations provide information about the evolution of isocurvature modes, which are 
controlled by the masses of the fields M as well as possible turns of the inflationary trajectory. 

The solution to the set of equations (2.6) or (2.11) can be formally expressed without 
any loss of generality in terms of a transfer matrix T(iV, W) which describes the evolution 
from time W at horizon exit to a subsequent time N 


dcf){N) 


= T{N,N,) 


S(f)* 


(2.13) 


“Note that after switching to the the orthonormal basis (2.9) there is no distinction between upper and 
lower indices. See Appendix A for a more detailed explanation of these expressions. 
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where the subscript (*) represents evaluation at horizon exit. Thus, all the information about 
the superhorizon evolution of the perturbations is encoded in the transfer matrix T. Using 
the decomposition Eq. (2.10), the transfer matrix takes a particularly simple form. Setting 
N = Nf 

=(ot) 1) (|i)' P''*) 

where the block T^_l has the same dimensions as a (n — 1) vector and T_l_l is an (n — 1) x 
(n — 1) matrix which represents the evolution of the entropy mode vector from horizon exit 
to the end of inflation. The matrix entry = 1 represents the requirement that curvature 
perturbations are conserved on superhorizon scales in the absence of entropy modes and 
T_l^ = 0 that curvature perturbations do not source entropy modes after horizon crossing [72] . 

In using the approach to computing ([, so far we have been taking the flat surface 

and constant density surface to be infinitesimally separated. A useful alternative is to take 
the flat surface to be at horizon crossing 


C(iVf) 


dN d4>i 


d(l){ d(j)% 

VtiV-d0 


0. 


(2.15) 


Written this way, all details of the superhorizon evolution are contained in the vector ViV. 
By direct comparison with Eq. (2.14) we can find and expression for ViV in terms of the 
components of the transfer matrix ^ 

= (2-16) 

We call the angle between ViV and the direction of gradient flow at horizon crossing ey the 
correlation angle A^. Since ey and T^_i_ are orthogonal to each other, using Eq. (2.16) we 
conclude that 

- e|| • eAT = cos Aat = (1 + (2-17) 

where T^_l = |T^_i_|, e^v is the unit vector in the direction of ViV, and the correlation angle is 
defined such that G [0, f ]• This quantity, as will become clear shortly, is a very convenient 
measure of the superhorizon evolution of C- 


2.3 Two-point statistics 

We define the power spectrum of the curvature perturbation to be 

27r2 

{C(k)C(k')) =(2^)=>i=>(k + k')^P(. (2.18) 

To make use of expression (2.15), we need to specify the conditions for the perturbations 
at horizon crossing. Provided the inflationary trajectory is not turning too much and slow- 
roll approximations hold so that H is approximately constant, it is reasonable to assume all 

^Despite the fact that the vector VA'^ contains information about the full super horizon evolution, it be¬ 
longs to the tangent space of M at the point (/)®(W), and therefore it is decomposed in the kinematic basis 
corresponding to the time of horizon crossing N = Nt. This is to be contrasted with Eq. (2.10), where the 
kinematic basis is that of = Nf. This implies a slight abuse of notation since we express the kinematic 
basis associated to these two different times with the same symbols Ba = {e||,e“}. 
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perturbations to be decoupled. Under these assumptions, the 2-point function for the field 
perturbations in a local frame where the fields have canonical kinetic terms, can be expressed 
as 


{d(^l{k)Scl>i{k')) = (27r)3ci*i53(k + k')^ , 


(2.19) 


where a star (*) indicates evaluation at the comoving scale k = aH. 

Applying the transfer matrix to the spectra at horizon crossing one immediately finds 
that the spectrum of curvature perturbations at the end of inflation has the form [62, 70, 71, 
73] 


^C = 




( 2 . 20 ) 


Since at horizon crossing T^_l = 0, according to Eq. (2.8), we see this expression can be 
simply written in terms of the correlation angle as 


Pc = 




COS^ 


A 


N 


( 2 . 21 ) 


Provided perturbations are small, this result relies solely on the assumption that the spectrum 
of perturbations is well described by Eq. (2.19) at horizon crossing. Note that an important 
consequence of assuming the field perturbations to be uncorrelated at horizon crossing is 
that to leading order in the the slow-roll parameters, superhorizon evolution of the scalar 
perturbations always gives a positive semidefinite contribution to the power spectrum. 

As we mentioned earlier, we can see from Eq. ( 2 . 11 ) that for cos Aw 7 ^ 1 we need sizeable 
entropy perturbations, i.e. isocurvature should not decay too fast, and we should also have 
mode mixing, that is, the effective mass matrix should have non-zero off diagonal terms. In 
the slow-roll regime the presence of mode mixing occurs whenever the inflationary trajectory 
deviates from the geodesic motion, or in other words, when it describes a turn in field space. 
In Eq. (2.11) these effects are encoded in the matrices M and Z. As we shall see in later 
sections these remarks have important consequences when trying to implement large scale 
suppression of the power spectrum in multifield inflationary models. 

We conclude this section by presenting an expression for the spectral index in this 
geometrical framework. As shown in Appendix B, in the slow-roll regime the spectral index 
can be written in a very compact way as 

dlogPr t - / ^ 

,, ^ =ns-l = -2e* + 2ej^ • M* • e^. (2.22) 


This expression is the same as that found in Ref. [58], which is the generalisation of the one 
presented in Ref. [62] for two-field inflation in the slow-roll slow-turn regime. Note that this 
result implies that the spectral tilt is determined by two local quantities, e and M, and one 
non-local quantity, the unit vector ew which depends on the details of the whole inflationary 
trajectory between the time of horizon-crossing iV* and the end of inflation Nj, since it 
depends on the transfer function T^_i_. When there is no superhorizon evolution the transfer 
matrix (2.14) reduces to the identity matrix, implying that T<^_l = 0 and ew is parallel to the 
direction of the inflationary trajectory ey. However, in general ew can point in any arbitrary 
direction with Aw determining the relative angle between ew and ey. 
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3 Power suppression in single-field infiation 


To understand the idea behind power suppression in the presence of a steepening of the 
potential, let us briefly discuss the case of single-field models. In this case, assuming slow- 
roll is a good approximation at horizon crossing, the power spectrum for the scale /c* = a*//* 
that crosses the horizon at is given by 


Pdk*) 



(3.1) 


Since during inflation the Hubble parameter is approximately constant, to realise P^{ki) < 
P(^{k 2 ), with ki < k 2 , reqnires that e(A^i) > e(A^ 2 )- In other words, power suppression 
between ki and k 2 can only occur if the slow-roll parameter e is decreasing between the times 
of horizon exit of these scales, Ni and N 2 respectively. This implies that the potential should 
be steeper between A^i and N 2 than in the subsequent inflationary evolution. 

The simple example proposed in Refs. [40, 45] illustrates this idea well. The potential 
has the general form 

V{^) = A^{Vs{cl)) + Vnm, (3.2) 

where Vs{4>) is the “slow” part of the potential, modelled as a Taylor expansion with coeffi¬ 
cients chosen to match observational constraints 


Hs = 1 - (3.3) 

The mass scale A is fixed by COBE normalisation [74] and makes Vs and lAi (and in particular 
(f)) dimensionless. The early stages of inflation are dominated by a brief period of “rapid” 
evolution on a steeper potential Vr, that shonld be steep enough to provide the desired power 
suppression. According to Refs. [40, 45], these requirements are satisfied by a quadratic 
potential like 

Vr = 0(())c - (/))^mi((/>c - (/>)^ (3.4) 

where (t)c can be chosen to produce power suppression at ^ < 50, as we do for example in the 
right plot of Fig. 1. 

The key observation made in Ref. [40] is that for the model given by Eq. (3.2) one can 
obtain an approximate expression for this power spectrum, when m\ ^ 1, of the form 

Pc«^|l-2m?^| Pi (3.5) 

where P^ is the power spectrum that would be obtained by considering only. Importantly 
this expression shows that the effect of lAi is always to suppress the power spectrum. This is 
intuitive, as provided Eq. (3.1) holds, one sees that the effect of the “rapid” phase manifests 
as increasing the magnitude of e. 

To understand more accurately the power suppression induced by the presence of Hr, 
one should take into account the deviations from the slow-roll approximation or from (2.19) 
induced by the period of rapid evolution. For this we performed a full numerical evolution 
of the background and two-point correlation functions of the flnctnations (which does not 
rely on slow-roll approximations) [55] For simplicity, and consistency with Refs. [40, 45], 

®Our numerical approach is the non-slow-roll version of the Transport method. We refer the reader to 
Refs. [55, 61, 75-77] for more information on this approach. 
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Figure 1. Power spectrum of the model (3.2) with parameters es = 0.005, A"* = 3.5 x 10“® and 
different values of the mass parameter m\. Both plots correspond to = 0.07, 0.1, 0.15, 0.3 repre¬ 
sented by a solid, long-dashed, short-dashed and dotted lines respectively. LEFT: We have set ko to 
be the scale k that exits the horizon at the time the Ng = 0. RIGHT: In order to make contact with 
observations it is more convenient to fix fco so that the scale k that exits the horizon at the time the 
transition occurs (j){Nc = 4) = satisfies \ogk/ko « 4. 


curvature effects were completely ignored in the evolution equations. We believe that given 
the current constraints on spatial curvature this approximation should not impact on our 
main results, which refer to scales leaving the horizon after the curvature bound has been 
reached. In addition, if our observable universe tunnelled from a parent metastable vacuum, 
and inflation started earlier than the horizon-crossing time for the largest observable scale 
in order to dilute curvature artefacts, it seems reasonable to assume a Bunch-Davies state 
as our initial conditions at this time. This is the starting point of our numerical evolution 
which we define to be at A^e = 0 The results of the simulations are displayed in Figs. 1 
and 2. 

In Fig. I we have displayed how the power spectrum varies for different choices of the 
mass parameter which determines the steepening, while the other two parameters in the 
potential, A and e, are kept fixed for simplicity. The spectrum in these plots interpolate 
between the spectrum of quadratic inflation for large scales k ^ ko and the one of linear 
inflation for small scales. The transition point corresponds roughly to the point where the 
spectral tilt — 1 changes abruptly from positive to negative. The figure shows that, as 
the “rapid” region of the potential becomes steeper (for increasing values of mf), the power 
suppression in the 4 e-folds prior to the transition becomes more pronounced, and at the 
same time the range of scales which experience power suppression decreases. 

In order to identify which models are capable of producing power suppression on large 
scales, it is useful to observe that the suppressed spectra necessarily have a positive (blue) 
tilt i.e. n* — 1 > 0, at some point during the first 4 e-folds of the observable inflation. We 
therefore use the deviations of the spectral index on large scales as a marker for the existence 
of power suppression. As we discussed in the previous section, the general expression for the 
spectral index in the multifield case is complicated by the fact that it involves non-local terms. 
However, in the case of single-field inflation, things become much simpler since the vector 

®The initial state after tunnelling is known to be affected by the parent vacuum as well as other effects 
associated to the nucleation process. Our assumption is that these effects are small for the cases studied in 
this paper. This point has recently been investigated in [49]. Their conclusions indicate that this would be a 
small effect. 
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Figure 2. Power spectrum of the model (3.2) with parameters es = 0.005 and = 3.5 x 10“®. We 
have set ko so that the scale k that exits the horizon at the time the transition occurs (j>{Nc = 4) = 
satisfies \ogk/ko « 4. LEFT: Power spectrum for the mass parameter = 0.07. The grey solid line 
represents the spectrum obtained from the slow-roll approximation (3.1) and the red solid line results 
from a full numerical simulation. RIGHT: Power spectrum obtained from a full numerical simulation 
with mass parameters mf = 0.07 (solid line) and mf = 2.4 (dashed line). 


e^v necessarily coincides with the direction of the inflaton ey, and therefore the expression 
for the spectral index to leading order in the slow-roll parameters, Eq. (2.22), reduces to the 
usual formula 

eJi • M • ey = —2e -|- r] — 1 = —6e -|- 2r]. (3-6) 

Here rj is the second slow-roll parameter of single-field inflation given hy t] = Using 

this expression for the spectral index we can translate the blue tilt condition on the spectrum 
into a constraint on the inflationary potential. According to Eq. (3.6), a blue tilt for the scale 
/c* occurs whenever 

r?>3e ^ v^^V>^-Vl. (3.7) 

We expect the slow-roll approximation to be in good agreement with the numerical result 
for scales exiting the horizon well away from the transition between the rapid and slow parts 
of the potential. However scales leaving the horizon during the transition between these two 
periods of inflation can experience effects that cannot be captured by the slow-roll description. 
One effect is the mixing of mode functions, which gives rise to oscillations in P(^{k) as can 
be seen in Eig. 2. This effect is well studied and we refer the reader to Refs. [36, 47, 48] for 
more detailed discussion 

In the left plot of Eig. 2 we compare the numerically obtained non-slow-roll result with 
the slow-roll approximation which we obtained solving the gradient flow equation Eq. (2.5) 
and assuming slow-roll at horizon crossing, Eq. (3.1). When the mass mi is small, as for 
the model parameters discussed in Refs. [40, 45], then the non-slow-roll correction is small. 
However, it is noteworthy that increasing the mass mi does not simply reduce the range of 
scales over which power suppression occurs; we also see that the oscillations became large. An 
example of this is shown on the right plot of Eig. 2, where we compare the spectra obtained 
with mf = 0.07 and mf = 2.4. If the potential is too steep the resulting power spectrum 

Recently, in Ref. [48], a study of the implications of different equations of state at the beginning of 
inflation was performed (for related work see [78, 79]). In this paper [48], authors considered an instantaneous 
change between pre- and inflationary eras. We believe that slow approaches to the inflationary attractor might 
also have interesting phenomenological consequences. 
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could be ruled out by observations, imposing important constraints on this model as a viable 
method to obtain power suppression. 


4 Power suppression in multifield infiation 

The construction of multifield inflationary models which display suppression of the power 
spectrum at large scales is more subtle than in single-field models. Consider the ratio of the 
power spectrum P(^{k) at two different scales ki < k 2 which exit the horizon at the times Ni 
and N 2 respectively. In order for the power spectrum to be suppressed at the scale ki with 
respect to /c 2 we need to have 

^ f cos^ Ajv, \ 

Pdk2) P^{k2) Vcos^A^J ^ ^ ^ 

In other words, this ratio depends on the relative amount of superhorizon evolution experi¬ 
enced by scales ki and k 2 , which in turn depends on the relative size of the corresponding 
correlation angles Aat^. Recall that Aat^ G [0 ,7r/2] is the angle between the gradient of the 
number of e-folds, parallel to e^r, and the inflaton direction at the time A*, ey. For a given 
ratio P^{ki)/P^{k 2 ) at horizon exit, we can find two different sitnations: 


• Aatj < • When the correlation angle is smaller for the smaller scale the superhori¬ 

zon evolution acts as to reduce the power suppression existing at horizon exit 


- pi{k2y 


(4.2) 


As a consequence of this, to have power suppression at the end of inflation in this 
scenario, it is not sufficient to have it present at horizon crossing. 


• Ajv* > Atv* : When the correlation angle is smaller for the larger scale the superhorizon 
evolution acts as to enhance the power suppression existing at horizon exit 


pyk2) ^ Pi{k2y 


(4.3) 


In this situation power suppression might be present at the end of inflation even when 
it is absent at horizon crossing, i.e. P^(/ci) « P^{k 2 )- 


For an inflationary model to be predictive, inflation must reach the adiabatic limit before 
the end of inflation. In scenarios where this is the case, it seems reasonable to expect that 
scales exiting the horizon at later times experience progressively less superhorizon evolution. 
From this perspective, the first situation seems somehow more typical than the second one. 
The second situation can occur, for example, when a very massive field suddenly becomes 
light at some point between Ni and A 2 . In that case, low power at large scales arises because 
of a power enhancement at small scales. In the next section we will show a concrete example 
of the first situation, and in § 6 we will present an example of the second. 
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Analysis of the spectral index: Similarly to the case of single-field inflation, a useful 
way to translate the conditions for power suppression into constraints on the form of the 
potential is to look at the evolution of the spectral index. As we argued above, the presence 
of suppression at large scales requires a period where the tilt transitions to blue. In the 
multifield case a blue tilt for the scale A:* occurs whenever 

ejy • M* • Cat > e*, (4.4) 

in other words, whenever the projection of the effective mass matrix M* along the vector 
ctv is larger than the slow-roll parameter e*. In general the vector e^r will not point in the 
direction of the inflationary trajectory and therefore, in contrast to the single-field case, the 
condition above is unrelated to the steepness of the potential. In fact, the steepening of the 
potential is neither necessary nor sufficient to guarantee a period of blue tilt in the large scale 
power spectrum. 

To understand the constraints on the potential around the inflationary trajectory, it is 
convenient to express Eq. (4.4) in terms of the Hessian of the scalar potential, rather than 
M which is related to the Hessian of log V: 

ejy • H* • otv > (1 + 2 cos^ AatJ K e*, where ffj = V^VjV - ^Rj. (4.5) 

O 

When the field space metric is flat everywhere the curvature term in the definition of H is 
absent and it reduces to the Hessian of V. Since the diagonal elements of any hermitian 
matrix are always bounded above by the largest eigenvalue, this condition implies that the 
largest eigenvalue of H* must be larger than at least Ke*. Therefore, a necessary but not 
sufficient condition to have a blue tilt in the spectrum at a given scale fc* is that the largest 
mass of the system evaluated at horizon exit satisfies 

^!i,|max ^ I4i, e*. (^-b) 

In the same way, since diagonal elements of a hermitian matrix are bounded below by the 
minimum eigenvalue, we find the sufficient but not necessary condition for blue tilt at 

^*|min ^ 3I4i,e*. (4-^) 

In other words, when this condition is satisfied the power spectrum will always be blue-tilted 
at k*, irrespective of the superhorizon evolution experienced by the mode after horizon exit. 


5 Example I: plateau preceded by quadratic inflation 


In this example we consider the situation where Atvj < ^Ni- We discuss a two-field model 
where the fields have canonical kinetic terms Gij = 5ij, and we define the scalar potential to 


be 


V{cPi, 02) = Hs( 0 l) + Eri( 0 i) + Hr2(02) 


(5.1) 


Here Vs and Hri remain the same as in Eqs. (3.3) and (3.4), and Hr 2 is given by 


kR2 = 



(5.2) 


This illustrates the picture where a mass hierarchy is present in the scalar sector resulting 
in non-trivial dynamics before the inflationary plateau. This is the typical situation after a 
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tunnelling event where the initial conditions for inflation are generically far away from the 
inflationary region of the potential. 

The first thing we can tell about this model is that in order to have power suppression, 
according to the necessary condition (4.6), we require that during some period 

max{mi,mi} > (5.3) 

where we have also used that in the steep region of the potential V > V((/>c) and e > es- 
If this condition is not satisfied, no slow-roll inflationary evolution will lead to a blue-tilted 
power spectrum, and therefore we can not have power suppression on large scales. 
Conveniently, this potential belongs to the sum-separable class 

Vict>i,...,cj>^) = Y,Viicj)i) (5.4) 


for which it is possible to calculate analytic expressions for observable quantities, provided 
the background evolution is well approximated by the slow-roll equations and the kinetic 
terms are canonical. Specifically, the vector VA^ in these models can be written as [80] 


(VfV)i|7V, 






,2 * 


(5.5) 


where Zi are a set of functions that become constant in the adiabatic limit. Without loss 
of generality we can always rearrange the constant terms of the summation V = Vi such 
that Zi\N^ = 0 Thus, using Eq. (2.20), we have that 


^C = 



(5.6) 


which depends exclusively on quantities evaluated at horizon crossing, even though it encodes 
all superhorizon evolution of the power spectrum. The term in curly brackets can be thought 
of as the result of transferring power from isocurvature perturbations to adiabatic perturba¬ 
tions on superhorizon scales. Using this expression and the condition for power suppression 
Eq. (4.1), we can immediately find constraints on the potential. The requirement for power 
suppression is not a period of inflation with relatively larger e (like in single-field inflation), 
but a period where the summation Yli (^/^*)^ is smaller^^, which is a significantly more 
intricate condition. 

To better understand how this translates to constraints on initial conditions it is con¬ 
venient to express Eq.(5.6) in terms of the correlation angle. Remembering that A^r is the 
angle between VU and ViV 


VU = A^(U5,i + Vri,i , mi <P2), VN = (5.7) 

a field reaches the minimum of its own potential during inflation = 0, it is always possible to 

redefine the potential such that =0 also vanishes at that point. It can be shown that this choice implies 

that the corresponding Zi has to be zero at (pg (see [80, 81]). 

^^N.B the explicit e factor in Eq. (5.6) cancels with the term in the denominator of . 


- 14 - 







Figure 3. Scenario One (mi = 0). We show several different inflationary trajectories on the (j)! — 4^2 
plane. None of the trajectories have power suppression at large scales even though it is present at 
horizon crossing. 


and that |Vy| = \/2eV^ we obtain 


cos A yr = 


1 

V^\VN\ 


< 

\h\ 


(5.8) 


When the system is in its adiabatic limit the two directions coincide and Ayr = 0. Conversely, 
for large values of 4)2 the correlation angle is close to tt/2, implying that in this regime the 
perturbations will experience a large amount of superhorizon evolution. Using the previous 
result it is straightforward to derive an upper bound for the magnitude of power suppression 
in the multifield case 


1 > 


Pdk2) 


> 


e\4)2 


Ni 


cos 


A 


N2- 


(5.9) 


For an analysis of initial conditions, it is convenient to make a change of variables in field 
space such that 

(pi = 4>c + ^ + R cos 9, (f >2 = Rsm9. (5.10) 

mf 


and start the trajectories always at the same height, from a fixed value of the scalar potential 
= bo- These initial conditions define a curve in field space which is characterised by a 
function R = R(9). Expressing 4)2 and e in terms of Vb and 0, and keeping only the leading 
order terms in the slow-roll parameter e, we can rewrite Eq. (5.9) as 


^ AUq eN2 cod An2 
~ (mf cot^ 9ni + m2) ’ 


(5.11) 


where AVb = Vb — is the initial height above the inflationary plateau. As we require our 
model to approach the adiabatic limit in order to be predictive, for simplicity we will assume 
that this limit is reached shortly after the steep region of the potential, such that Ayrj = 0 . 
Moreover, we will also suppose that at this point the inflaton is rolling down the inflationary 
plateau where = 65 . In this case, the condition for power suppression can be expressed 


- 15 - 














Figure 4. Superhorizon evolution of the power spectrum for scenario one — the two-field model 
given by Eq. (5.1) where = 0 and m| = 0.2. LEFT: P({k) changing from horizon exit until the 
end of inflation. RIGHT: comparison between P(^ evaluated at horizon crossing (dashed grey) and at 
the end of inflation; lighter red is the numerical result and dark red is the analytic expression given 
by Eq. (5.6). We see that initially, at iV = 0, there is signihcant power suppression on large scales 
(small \og{k/kQ)) but soon after the largest scale leave the horizon, the inflationary trajectory goes 
through a turn, thereby giving rise to superhorizon evolution which acts to erase any information of 
this initial state, fcp, A and es are defined as in Fig. 2. 


in terms of a maximum value of the initial angular direction, 0|max = dc, above which power 
suppression cannot be realised: 

2 /I 

tan^ Oc = ——-7^-^ 

es AVo — A^m2 

Scenario One (mi = 0) : The effects of isocurvature are especially evident in the case 
mi = 0. As argued in § 4, in order to have power suppression at horizon crossing the 
inflationary trajectory needs to experience a steepening, and in this case the rapid evolution 
occurs in a different field direction to the slow phase, i.e. along (j) 2 - This naively looks like 
a small change from the single-field scenario since one can still arrange for the trajectory to 
follow very similar values of the potential thus resulting in a similar profile for e. However the 
role of multifield effects turns out to be very important - when the superhorizon evolution 
is taken into account the power suppression disappears, completely. This follows directly 
from Eq. (5.12), as the critical angle vanishes and therefore there are no initial conditions for 
which one can realise power suppression. 

The evolution of the power spectrum on each scale from horizon exit up to the end of 
inflation is shown in Fig. 4. Around horizon crossing, one sees significant power suppression 
on large scales (for logk/ko < 3) but soon after these scales leave the horizon, the field-space 
trajectory goes through a turn (see Fig. 3), causing the isocurvature perturbation to source 
the adiabatic perturbation and hence evolves in such a way as to erase any information 
about the steepening of the potential that was there initially. 

The change in the spectral tilt displayed in the right plot of Fig. 4 can be understood 
easily from the multifield expression of the spectral index Eq. (2.22). Consider a trajectory 
starting at some large value of (f >2 — initially ey points roughly along the 4>2 direction and 
e_L points along the (/)i direction (see Fig. (3)). If there was no superhorizon evolution, the 
spectral index would then reduce to 

n* - 1 -6e* H-——, (5.13) 


(5.12) 
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Figure 5. Scenario Two {m 2 , mi > 0). The parameters of the model are set to es = 0.005, = 

3.5 X 10“® and ml = 0.07. The second mass parameter is given, from top to bottom, by = 

ml, 0.25 X ml, 4 x ml- LEFT: Inflationary trajectories on the (j)i — (j )2 plane. The trajectories 
associated to a suppressed power spectrum are drawn in red. RIGHT: Power spectra corresponding 
to the initial conditions represented by blue dots in the left figures. From bottom to top the solid, 
long-dashed, short-dashed and dotted lines correspond to increasing initial values for (j )2 at N = 0. 
To see more clearly the effect of changing the initial conditions, ko is fixed as in the left plot of Fig. 1. 


implying that at a certain height of the potential V we could have a blue-tilted spectrum 
by choosing to be sufficiently large. This would be the naive conclusion one would get 
by studying the spectrum at horizon crossing. However, making superhorizon evolution into 
account, from the bound (5.8), we see that at large 1^2 the correlation angle should be close 
to 7r/2, implying that e^r points instead along the direction of (/>i. As a consequence, using 
Eq. (2.22) with = 0, we find that the spectral index is smaller than one 


Us-l^ -2e*, (5.14) 

i.e. the spectrum is red-tilted at the end of inflation. This comparison is shown explicitly in 
the right plot of Fig. 4. 
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Figure 6. Scenario three with mf = 0.2 and = 40 x 0.2. LEFT: field trajectory; RIGHT: 
corresponding P({k). The power suppression closely resembles the single-field example, except now 
there are also smaller scale oscillations resulting from the oscillations of the second field (j) 2 - kg, A 
and es are defined as in Fig. 2. 


Scenario Two {m 2 ,rni > 0) : In our second scenario we consider the effects of having 
a non-zero value for mi. Taking mi to be the same as in the left-hand plot of Fig. 1, i.e. 
appropriately chosen to manifest power suppression when the second field is placed at its 
minimum, we conclude that the results depend heavily on the choice of initial conditions. 
In Fig. 5 we show the inflationary trajectories and corresponding power spectra of the cases 
mi = m 2 , mi > m 2 and mi < m 2 respectively. In these cases depending on the choice of 
initial conditions the result interpolates between the full suppression of the single-field model 
and no suppression, when ( 1)2 plays a significant role in the evolution, as occurs in the previous 
example. 

The left-hand plots show the slow-roll inflationary trajectories for each case. The red 
lines correspond to inflationary trajectories where the corresponding power spectrum displays 
suppression. For convention, we say that a power spectrum has suppression when the deficit 
of logP^ four e—folds before transition with respect to the transition itself N = Nc 

logP^lAT^ - logP^|Ar,-4, (5.15) 


is at least half of what it is obtained in the single-field case. The plots illustrate that the 
angle Oc increases with the magnitude of m 2 , as expected from the formula (5.12). 

Scenario Three (m 2 mi) : For completeness, our third scenario is the same as the last 
one, except now m 2 mi. This case cannot be treated using the slow-roll approximations 
since displacement of the second field can give rise to significant oscillations. As expected 
from basic considerations of effective field theory, as the limit m 2 S> mi is approached, the 
situation becomes simpler and the broad features match that of the single-field result. The 
difference is that now we also see oscillations on the largest scales as a result of excitation 
of the second field; see Fig. 6. This last scenario displays essentially the heavy physics 
phenomenology that has been studied in, for example. Refs. [82-88] and references therein. 

The nature of the oscillations of a second field are different to that of the mixing of 
modes observed in the single-field scenario. It was shown in Refs. [84, 85] that the features 
related to the decay of a heavy mode behave as 


ZAP 


osc 


oc —I COS 



k 

Cl:^7Tih 


(5.16) 
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Figure 7. Scalar potential in flume inflation, defined by equations (6.1) and (6.2). The inflaton rolls 
down the potential mainly along the angular direction. The mass associated to the radial direction, 
m^, changes abruptly from a heavy phase where ^ H^, and thus the model is effectively single 
field, to a light phase where < H^, and therefore the effects of superhorizon evolution of the 
perturbations become relevant. 


where ruh is the mass of the heavy field and Pq is the unperturbed power spectrum. The 
mixing of modes as a result of a sudden change in the background evolution [36, 89, 90] gives 
rise to a different oscillatory profile. Even though the amplitude of oscillations is sensitive 
to details of the period of rapid evolution, in the short wavelength limit it is possible to 
characterise their frequency and decay rate in a relatively model independent way. 

Provided the density perturbations asymptote to the Bunch-Davis vacuum for scales 
deep inside the horizon, the power spectrum has the following characteristic behaviour for 
k/ttcH ;§> 1 (see Ref. [89]) 



(5.17) 


where Oc is the value of the scale factor at the transition N = N^- With sufficiently high 
quality data one could hope to distinguish these two signatures. 


6 Example II: Power suppression from superhorizon evolution 

So far all our examples have been simple extensions of a single-field model with an initially 
steeper phase of inflation. In the examples shown, multifield effects only act as to reduce the 
amount of power suppression since larger scales experience more superhorizon evolution than 
smaller scales. However, as described in § 4, in principle it should be possible to generate 
power suppression on large scales purely as a consequence of smaller scales experiencing a 
greater level of superhorizon evolution. This is a very different mechanism to that of the 
previous sections and, to our knowledge, there are no examples of models of this kind in the 
literature. The purpose of this short section therefore is to provide such an example. 

The model we propose is inspired by inflation in not single-valued potentials which are 
present in descriptions of moduli spaces in string theory models of compactifications (see 
also [91, 92]). However, for the purposes of this paper, we will treat it purely as a toy model 
and make no claims of a particle physics embedding. This interesting possibility is left to 
future work. To be getting on with, we shall simply refer to this preliminary model as Flume 
Inflation. The functional form of this model is very similar to the previous examples, except 
now we work in polar coordinates {r, 9}. The field-space metric is thus 

Gii = (J " 2 ) , (6.1) 
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Figure 8. The lefthand plot shows the trajectory to spiral outwards as result of the decreasing radial 
mass. The righthand plot demonstrates that while the spectrum is close to flat at horizon crossing, 
superhorizon evolution causes the amplitude on smaller scales to grow signiflcantly, resulting in a 
spectrum which is more commonly referred to as “suppressed” on the largest scales. 


but note that the field space JW has zero curvature, and therefore geodesics are still straight 
lines in cartesian coordinates. The potential has the form 


V = A^ 


1 + + ^m(0f(r - ref , 


( 6 . 2 ) 


plotted in Fig. 7, where the key step is to make the radial mass a function of 0, such that 
it becomes lighter at some point 0c while relevant scales are leaving the horizon. To achieve 
this in an intuitive way we use the error function 

m(0f = + + erf ) (^i “ ^h) • (6-3) 

where determines the length scale over which the mass changes from m to m;. Initially we 
set ~ such that isocurvature modes exiting the horizon during this phase are 

rapidly decaying and hence the evolution, despite the curved trajectory, is essentially single 
field. After the transition, the mass rapidly approaches mf <C such that the slow-roll 
approximations are valid at late times and hence the model is amenable to the discussion in 
§ 4. During the transition from the heavy phase to the light phase the approximations of 
§ 4 do not hold and hence the results shown are obtained numerically without making any 
slow-roll approximations. Nevertheless at a qualitative level we will see that the intuition 
obtained from § 4 serves us well. 

The dynamics of the inflaton is reminiscent of a person riding down a spiralling flume. 
Much like the person, the inflaton will sit close to the minimum in the radial direction, 
and will roll mainly along the angular direction. This deviation from the geodesics causes 
excitation of isocurvature modes and also transfers power from isocurvature perturbations 
to adiabatic perturbations, but this effect is only significant when the radial mass becomes 
sufficiently light. We therefore expect Ajv 7 ^ 0 for scales exiting the horizon during the later 
stages when m? ~ mf. Hence, by arranging for the radial mass to become lighter during 
the period when CMB scales are exiting the horizon, we ensure that smaller scales undergo 
a larger amount of superhorizon evolution than larger scales 

this behaviour continues until the end of inflation, then the adiabatic limit will not be reached and 
hence the model is in principle sensitive to the details of reheating. We will not address this possibility here. 
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An example of this is shown in Fig. 8. The lefthand plot shows that although the 
minimum is at a fixed radius, the trajectory spirals outwards as a result of the decreasing 
mass. The righthand plot shows that at horizon exit the power spectrum is close to flat. This 
is a direct consequence of e ~ constant. We then see that smaller scales undergo a substantial 
amount of superhorizon evolution while large scales do not. This is because the isocurvature 
on the largest scales underwent a period of exponential decay and hence, even when the radial 
mass becomes light, isocurvature modes never grow to the extent that isocurvature modes 
exiting the horizon during the lighter phase do. The result is that by the end of inflation the 
spectrum does indeed appear suppressed on large scales. 


7 Discussion 

Observing statistically significant power suppression on the largest observable scales would 
require some non-minimal version of inflation and therefore radically improve the possibility 
of extracting valuable information about fundamental physics from the CMB. A number of 
authors have pointed out that a steepening of the potential associated with a tunnelling 
event can result in suppression of the power spectrum on large scales. Our paper extends 
this analysis by discussing the implications of this in the context of multifield inflation. 

We find that the existence of a steepening along the inflationary direction is not a 
sufficient condition for power suppression, but rather that there is a relevant direction in 
field space that determines whether or not the power spectrum presents this suppression. 
This direction depends on the full inflationary trajectory in field space and is therefore not 
something determined locally from the potential. Furthermore, its orientation does not have 
to be aligned with the velocity vector in field space so fast evolution does not necessarily lead 
to power suppression even if it is present at horizon crossing. The reason for this behaviour 
is that the superhorizon evolution of isocurvature modes can source the adiabatic mode, 
thereby restoring its level of power after horizon crossing. 

For models where the potential can be expressed as sum-separable, the situation becomes 
sufficiently simple such that it is possible to give analytic expressions for the constraints on 
the potential at horizon crossing to ensure large scale suppression of the power spectrum. This 
is further simplified in models where the adiabatic limit is reached shortly after the transition 
from the steep to the flat part of the potential where most of the e-folds of inflation would 
take place. We present simple 2-field models where we can easily visualise the dynamics and 
show explicit examples where the suppression at horizon crossing completely disappears after 
superhorizon evolution is taken into account. We also study the resultant power spectrum 
as a function of the parameters of the potential and initial conditions. This is important in 
the landscape picture where the initial conditions would be determined by an instanton that 
interpolates between the previous vacuum and the one described by this potential. 

We also describe a novel approach to power suppression where the superhorizon evolu¬ 
tion of perturbations plays a key role. One can envision a model where the relative power 
suppression is created by an enhancement of the level of the perturbations after the first 
few e-folds of inflation, in other words, one does not suppress the power at large scales but 
enhances it at small scales. We present a toy model that achieves this as a result of transfer 
of power from isocurvature to adiabatic perturbations over a prolonged period. This exam¬ 
ple may have interesting applications in inflationary models arising from compactification 
scenarios. 
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Our results have interesting implications for inflation in the landscape. If our observable 
patch emerged after tunnelling from a metastable vacua, then the combined conditions on 
the potential of the active helds and their initial conditions in order to achieve a suppressed 
large scale power spectrum imposes strong constraints on where in the landscape this event 
might have happened. In summary, observation of power suppression on large scales would 
enable one to establish a connection between the initial conditions for the onset of inflation 
in the landscape and the dynamics of active fields. We have made first steps in making this 
relationship precise. 
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A Evolution of the scalar perturbations in the kinematical basis 

In this appendix we will derive the equation for the scalar perturbations on superhorizon 
scales in the kinematic basis, where the evolution of curvature and isocurvature perturbations 
is explicit. We follow the notation of Refs. [62-64] . In these works, the curvature perturbation 
is defined in the comoving gauge and therefore denoted TZ. Curvature perturbations dehned 
in the constant density gauge — Q — and in the comoving gauge — TZ — are known to 
be equal at second order on superhorizon scales up to 0{k/aH)‘^ corrections [67, 68]. Eor 
consistency with the standard usage of the separate universe assumption, in this paper, we 
will always refer to the curvature perturbation as C,. 

Eor simplicity we will just discuss the case where the held space metric is hat everywhere, 
and therefore will not make any distinction between upper and lower indices. Moreover, we 
restrict ourselves to the slow-roll regime, where the background satishes the equation 

0' = -VlogU. (A.l) 

In particular we will show that curvature perturbations are conserved on superhorizon scales 
in the absence of entropy modes, and that curvature perturbations are not transferred into 
entropy modes. This two facts were used in § 2 to explain the simple structure of the transfer 
matrix T in (2.14). 

A.l Equations for the scalar perturbations 

As we discussed in § 2, scalar perturbations <50 on super horizon scales satisfy 

= -M ■ Scj) where M = VV log V. (A.2) 
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The evolution of scalar perturbations becomes particularly simple when the perturbations 
are decomposed in the kinematic basis, Bk = {e||,e"}, which is adapted to the background 

evolution. The basis vectors ey and (we use (2) to refer to the element 0 = 2, not to be 
confused by an exponent 2) are defined as 


where v = (A.3) 

= I </)'!, and Z 21 is defined to be a real 
positive number. The precise definitions of the rest of the basis vectors e“ with a > 2 
will not be necessary for our derivation, and can be found in Ref. [64]. For our purposes it 

suffices to know that the full set forms an orthonormal basis, i.e. = 5ab^ where 

i2'] 

Ga = {G||,e“}. As in the case of the vector ej_ , the evolution of the rest of the basis vectors 
along the inflationary trajectory is encoded in the elements of matrix Z 


eii = 




De\\ 

~dN 


= Z2ie 


( 2 ) 


The vector eii has unit norm since v = |VlogR| 


^ab — 6a 


-Pefc 

dN ■ 


(A.4) 


It is easy to check that this matrix is antisymmetric from the condition ■ e^) = 0, 

which follows from the orthonormality of this set of vectors. In the kinematic basis the 
perturbations Sip can be decomposed as 


N 

S(f) = uCeii + ^ (5(/>" e“. 
0=2 


(A.5) 


As we mentioned above, equations (A.2) become particularly clear when the perturbations 
are decomposed in this way, and this is due to the simple structure that the matrix M has 
on the basis Bk- To see this first note that from the background equation (A.l) we can derive 
the following two relations 

^ = 0'.V = -(VlogR)-V ^ ^ = -^iy\ogV) = -M-<P'. (A.6) 

Then, using cp' = uey and the definition for the basis vector (A.3), we can rewrite the 
second equation as 

Projecting this relation along the vectors = {e||,e"} we find that the elements of the 
matrix M in this basis satisfy 

eij • M • eii = —• M • en = —Z 21 , • M • en = 0 for all j3 > 2. (A.8) 

11 II y _L II _L II 

Finally, when the expression for the perturbations in the kinematic basis (A.5) is substituted 
into equation (A.2) we find 


dv 

dN 


Cey +v 


dc 

-e 

dN 


+ vC 


De\\ 

'm 


+ 


dSpl 

dN 


+ 6PI 


Del 


-vCM • ey -SplM ■ e([. 


(A.9) 
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and after projecting this expression along the vectors ey and e“, we hnd that the curvature 
perturbations ( and isocurvature perturbations (/>" satisfy respectively 


dN 

dN 


—2Z21 



V 




(A.IO) 


From this equations it is now clear that curvature perturbations are conserved in the absence 
of isocurvature pertubations, and that isocurvature is not sourced by curvature perturbations. 
It is straight forward to check that, when the second equation is written in terms of the 
entropy perturbations 5" = it takes the form 


In 




where Muh = en • M ■ e 


(A.ll) 


A.2 The transfer matrix 

The set of equations (A.IO) and (A.ll) can be written in matrix form as follows 


dAT \Sl‘) ^ ■ ' 

The transfer matrix T{N 2 , Ni) between two times labelled by the e-fold numbers Ni < N 2 
is defined as 

(A,13) 

And therefore, the equation for the perturbations (A. 12) implies that the transfer matrix 
between the times N, and N + dN is given by 


T{N + dN, N) 


(1 — 2.^21 < 52/3 

\^0 (1 - M||||) - [M - Z]ap 


Since, by definition, transfer matrices satisfy the composition rule 


(A.14) 


T(A3, Ai) = T(A3, A 2 ) • T{N2, Ai), (A.15) 

the full transfer matrix T{N 2 , Ni) can be obtained formally by multiplying a sequence of 
infinitesimal transfer matrices of the form (A.14). Moreover, as the composition preserves 
the upper triangular structure in (A.14) and also that T<;/^ = 1, it follows that the structure 
of the full transfer mass matrix T{N 2 , Ni), has to be of the form 

T=(jjy). (A,16) 

B Spectral index in multifield inflation 

For convenience we define the transfer matrix T as 

f=—T dct){N) = f{N,N^)-Sct,^, (B.l) 

V 
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which satisfies the same formal equation as Scf) 


DT 


{N,N,) 


-M f, 


(B.2) 


Transfer matrices T satisfy the composition rule (A. 15), which in particular implies that 


f{N,N^)-f{N^,N)=I. 


Taking a derivative with respect to N^, of the previous expression we find 


Df{N,N^) 

dN^ 


fiN„N) 


-f{N,N,) 


Df{N^,N) 

dN^ 


(B.3) 


(B.4) 


Then, using the equation for T (B.2), and projecting along the inflationary direction at the 
end of inflation, eidA^) we arrive to 


DT 


{N,N^)=f-M, 


DTi 




dN* ^ ’ dN* 

where we have defined T|| = e||f(A^) • T. Finally, writing the power spectrum as 


C = -Til • Sc/)^ 


P( = 




(B.5) 


(B.6) 


and given that the first equation also implies ejv = T|/|T||,we have 


d log Pq 
dN* 


— —2e* + 


2 DP\\ 

Ti I dN^ 


gtv — —2e* + • AT • gtv- 


(B.7) 
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